Theoretical evidence for efficient p-type doping of GaN using beryllium 
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Ab initio calculations predict that Be is a shallow ac- 
ceptor in GaN. Its thermal ionization energy is 0.06 eV in 
wurtzite GaN; the level is valence resonant in the zincblende 
phase. Be incorporation is severely limited by the formation 
of BesN2. We show however that co-incorporation with reac- 
tive species can enhance the solubility. H-assisted incorpora- 
tion should lead to high doping levels in MOCVD growth after 
post-growth annealing at about 850 K. Be-O co-incorporation 
produces high Be and O concentrations at MBE growth tem- 
peratures. 

PACS numbers : 71.25.Eq, 61.72.Vv, 61.72.Ss 



Gallium nitride is a base material for green-to-UV op- 
toelectronics, and piezoelectric and high-power devices. 
For these applications, controllable doping is an obvious 
necessity. Sic a is a suitable shallow donor, but the quest 
for the optimal acceptor still continues. The current-best 
p-dopant jl| is Mgc a , with which hole concentrations in 
the 10 17 cm -3 range have been achieved at room tem- 
perature (despite its ionization energy of ~ 0.2 eV), and 
whose activation mechanism is still under scrutiny [0. 

We show here that substantial improvements can be 
achieved by doping GaN with Be: ab initio calculations 
of formation energies, impurity levels, and atomic ge- 
ometries of Be in GaN show that Be is the shallowest 
acceptors reported so far in GaN. Calculated solubilities 
indicate barely acceptable doping levels in N-rich growth 
conditions at high growth temperatures, while complete 
compensation occurs in Ga-rich conditions due to N va- 
cancies. We show however that the co-incorporation of 
Be with H or O strongly enhances its solubility (other- 
wise limited by the formation of Be nitride) and may lead 
to p-type doping. 

Method - By means of density-functional-theory 
(DFT) H ultrasoft-pseudopotential plane-wave calcula- 
tions of total energies and forces in doped GaN (details 
are reported elsewhere ) , we predict from first prin- 

ciples the formation energies and thermal ionization en- 
ergies, and ensuing dopant and carrier concentrations of 
Be in GaN. The carriers concentration at temperature T 
due to a single acceptor with thermal ionization energy 
e[0/-] is 



D = N s exp (-£ form /k B T g 



(2) 



for a growth temperature of T g , with N s — 4.33xl0 22 
cm~ 3 available sites, and a formation energy -Eform- 
Omission of the entropic factor exp (Sj/ke) (formation 
entropies are currently impossible to calculate reliably) 
does not alter our conclusions. The formation energy for 
Beca in charge state Q is 



E toTm (Q) = E tot (Q) 



X 



+ Q{He + E$), (3) 



= D exp (-e[0/-]/k B T) = D£. 



(1) 



The impurity concentration D in thermal equilibrium is 



where fi e is the electron chemical potential, E tot (Q) is 
the total energy of the defected supercell in charge state 
Q, its top valence band energy, n x and /z x the num- 
ber of atoms of the involved species (X=Ga, N, Be) and 
their chemical potentials. The latter potentials must sat- 
isfy the equilibrium conditions with GaN and Be-X com- 
pounds. To obtain maximum solubility, we will always 
be choosing the highest /i Bc compatible with the rele- 
vant solubility limit. Then, a single chemical potential 
(e.g. /i N ) remains to be determined by the imposed (N- 
rich <-> Ga-rich) growth conditions. Since T=0 in our 
calculations, \x e is synonimous with the Fermi energy. 

The thermal ionization level e[0/— ] of a single acceptor 
is by definition the formation-energy difference of the two 
charge states Q=0 and Q=-l at /i e =0. Thermal ioniza- 
tion corresponds to hole release from (i.e., electron pro- 
motion into) the acceptor extrinsic state with concurrent 
geometric relaxation. Optical transition energies can be 
obtained as e.g. E gap — e[0/— ]— £"fc f° r conduction band- 
to-acceptor recombination. Epc is the Franck-Condon 
shift for the Q=-l state in the equilibrium geometry of 
the Q = state. For Be, E FC = 0.04 eV [|. 

The above quantities are DFT total energy differences 
and as such exact in principle, were it not for the lo- 
cal density approximation 0] and other technicalities. 
Among the latter, the one introducing the largest un- 
certainty is the alignment of the energy zeroes for the 
different charge states. The resulting error bar is typi- 
cally ±0.1 eV in the ionization energies, but in fact lower 
(~ 0.05 eV) for Be Ga (see Reds. §§ for details). 

Results - The thermal level of BeQ a is 0.06 eV in 
wurtzite GaN. This value agrees well with our previous 
estimate || based on a different (and less accurate) tech- 
nique. In the zincblende phase, the level is found to be 
resonant with the valence band. The difference between 
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zincblende and wurtzite is due to symmetry-related elec- 
tronic structure features, namely the impurity state being 
mostly localized on all the four impurity-neighboring N 
atoms in zincblende, while only on the a-plane neighbors 
in wurtzite (see Ref. M for more). Despite this minor 
deviation, the overall energetics is essentially the same in 
the two phases (the numbers below are for wurtzite). 

Beca is a shallow dopant in GaN, with an activation 
efficiency £ ~ 0.1 (sec Eq.[l]) at room temperature. To 
evaluate its chemical concentration, we calculated its for- 
mation energy as a function of growth conditions and of 
the Fermi energy (accounting for the dependence of the 
latter on charge neutrality conditions |(| involving back- 
ground doping and competing defects). For neutral Be, 
E imm is 2.29 eV (2.83 eV) in N-rich (Ga-rich) conditions; 
for fi e > e[0/— ], the negative state is favored (-Ef or m=2.35 
eV - fJ- e ). These high formation energies stem from the 
strict solubility limit imposed by the formation of Be 3 N 2 
(whose total energy is also calculated ab initio). There- 
fore, the Be concentration achievable in GaN is typically 
quite low. Also, competing defect configurations exist: 
we considered the nitrogen vacancy Vn, two kinds of in- 
terstitial Be (in the wurtzite channels and in the "trigo- 
nal" cage), and heteroantisite BeN (the latter being ruled 
out by a formation energy of ~ 8 eV) . 

The most favorable growth condition for BeQ a incor- 
poration is N-rich. is an ionized single donor with 
a formation energy of 2.32 + /i e (possibly behaving as 
multiple donor in the extreme p limit Q). The Be in- 
terstitials are double donors over most of the relevant 
Fermi level range (since e[2 + /+] ~ 2.3 eV), with for- 
mation energies of 2.39 eV + 2 \i e (cage) and 2.43 eV 
+ 2 [i e (channel). Therefore, none of these defects is 
ever competitive with Beg a (£f orm =2.29 eV) or Beg a 
(-£form=2.35 eV - fx e ). In particular, compensation by 
is suppressed. However, at typical MOCVD growth 
temperatures of 1500 K, the maximum concentrations 
reached is only D ~ 10 17 - 10 18 /cm 3 assuming formation 
entropies between zero and 5 ke (the entropic factor en- 
hances the solubility) . In wurtzite GaN at room temper- 
ature, thanks to Be's shallowness, this implies a doping 
level of 10 16 -10 17 cm~ 3 . For comparison, the calculated 
concentration H] of MgQ a in the same conditions is about 
10 19 /cm 3 (£form=l-40 eV), but the carrier concentration 
is low because of £ ~ 10 -4 (Eq.|l|) at room temperature, 
as implied by the calculated ionization energy of 0.23 eV. 

In Ga-rich conditions the situation is more critical. 
The formation energy of is 0.67 eV + /i e , whereas 
that of BeQ a is 2.89 - /x e . This results in complete com- 
pensation, with the Fermi level pinned at about 1.1 eV 
above the valence band top. The interstitials, again dou- 
ble donors with formation energies of 1.29 eV + 2 [i e 
(cage) and 1.33 eV + 2 /i e (channel), are competitive 
with BeQa only in p-type conditions (Fermi levels below 
~ 0.4 eV). Their action is anyway preempted by V^, 
which pins the Fermi level much higher in the gap. 



Co-incorporation - As shown earlier, Be alone could 
never be incorporated massively in GaN. An abundant 
supply of Be in the growth process would eventually cause 
Be3N2 to form, with detrimental effects for GaN crystal 
quality ffi] . A possible way to enhance the solubility of Be 
in GaN is "incorporation catalysis" by reactive species. 

We consider first the H-assisted incorporation sug- 
gested recently Q for Mg. The formation energy of a 
(neutral) Be-H complex in its equilibrium configuration 
in GaN is 0.32 (0.87) eV/complex in N-rich (Ga-rich) 
conditions (solubility limits: Be3N2 and H2 molecules): 
the concurrent incorporation of Be and H in GaN is 
therefore efficient. Accounting for competing species and 
charge neutrality, we find that as much as 10 19 to 10 21 
Be's and H's per cm 3 get incorporated in the Ga-rich 
and N-rich limits respectively, at a growth temperature 
of 1500 K. 

However, the acceptor is compensated by H, which 
must be eliminated from the sample to activate p-type 
conduction. After H is detached from the complex, Be is 
still compensated, since H is in its positive charge state 
||||. H should thus be removed from the p-doped re- 
gion by annealing in e.g. a N2 atmosphere. Being quite 
stable, the Be-H complex is easy to incorporate but dif- 
ficult to break up: the minimum binding energy differ- 
ence between Be-H and isolated BeQ a plus isolated H + 
in its equilibrium configuration (N backbonding site) is 
1.8 eV We estimate the barrier for detachment to 

be at least of order 2.1 eV. Within a simple desorption- 
like description, the annealing temperature T ann needed 
to break up H-Be complexes may be estimated from the 
rate 

R ~ i^CXp (--Ebarr/keTann). 

The relevant attempt frequency v is that of the H-Be 
bond center wagging mode, since the complex ground 
state is H at a Be-N bond center ||. We have previ- 
ously calculated || the stretching frequency for the H- 
Mg complex to be about 3000 cm -1 in the equilibrium 
N-backbonding site; as this Mg-H frequency essentially 
originates from N-H pairing [g|8j, it is expected to be 
similar for H-Be. As wagging frequencies of H at bond 
centers are usually about 10 times lower than stretching 
frequencies [0, we pick ^=100 THz - 300 cm" 1 . With 
these parameters, massive break-up of H-Be sets in at 
about 850 K, which seems a practicable post-processing 
temperature. 

Be-H co-incorporation seems possible. The analogous 
idea of reactive co-doping |ll[ also appears to be useful. 
In recent experiments [O], Be and O concentrations in 
the 5xl0 20 cm -3 range have been reported in zincblende 
GaN grown by MBE at 950 K in moderately N-rich 
and oxygen-contaminated conditions. Importantly, ex- 
periment jll| indicates that Be prevails over O by about 
20 % (D[Be]= 5 x 10 20 cm" 3 , D[0}= 4 x 10 20 cm" 3 ), 
leading to p-type doping. 
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We find that the co-incorporation of Be and O favors 
the precipitation of BeO over the formation of the indi- 
vidual solubility- limiting compounds BesN 2 and Ga2C>3, 
because the reaction 



Be 3 N 2 + Ga 2 3 -> 2 GaN + 3 BeO 



(4) 



is strongly exothermic (6.4 eV). Therefore, Be and O can 
be incorporated concurrently in large amounts into GaN. 
However, the incorporation of Be in excess ofO (as it ap- 
pears to occur in the co-doping experiments pJ|) is still 
affected by the BesN 2 solubility limit, so that the net 
Be concentration cannot be larger than that obtained in 
the absence of O. In consideration of the above-cited ex- 
perimental evidence by Brandt et al. |llf| , our results, 
obtained assuming the validity of equilibrium thermody- 
namics, imply that non-equilibrium effects are operative 
in the actual co-doping process. Such as-yet unknown ef- 
fects may be related to growth kinetics, as is quite plausi- 
ble considering the highly non-equilibrium nature of the 
MBE growth process. 

One may speculate that the presence of high chemical 
concentrations of Be (albeit O-compensated) will make 
it more likely for small non-equilibrium concentration 
fluctuations to occur in an interesting range (say ~10 19 
cm" 3 , i.e. ~ 5 % of the total concentration), whatever 
the exact cause of such fluctuations. This hypothesis is 
expecially intriguing, since the same mechanism may ap- 
ply to III-V nitrides in general: we found in fact that 
exothermic reactions analogous to Eq. (^) hold for sev- 
eral Group-II dopants (Be, Zn, Mg, Ca) in A1N and InN 
as well as in GaN, so that O-acceptor co-incorporation 
will occur efficiently in most cases. A treatment of this 
issue is beyond the scope of this work, and is deferred to 
a future work. 

We note that calculations show also that a neutral BeO 
pair substituting for a GaN pair is bound by 1.0 eV/pair 
with respect to over isolated Beg a and Oj: this path 
for O-Be donor-acceptor compensation may thus be pre- 
ferred over hole-electron recombination. BeO pairs carry 
at most a dipole, which scatters carriers less efficiently 
than a charged center. An enhanced carrier mobility 
should ensue, in qualitative agreement with unusually 
high mobilities being observed in experiment [ pd| . 

In summary, ab initio calculations predict that Be is 
a shallow acceptor in GaN. Be solubility in GaN is low 
due to Be3N 2 formation, so that hole concentrations at 
room temperature of only ~10 17 cm -3 are obtained from 
N-rich, high-temperature growth. In Ga-rich conditions, 
vacancy-acceptor compensation produces semi-insulating 
material. The presence of H enhances Be solubility in 
GaN, although acceptor activation requires annealing at 
about 850 K. Our finding of massive co-incorporation of 
O and Be in GaN grown in an O-contaminated atmo- 
sphere confirms in part the results of recent experiments 
Iffi!. In conclusion, further consideration of Be in actual 



growth experiments may open the way to a more effective 
p-doping of GaN. 
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